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We consider a dilute gas of neutral unpolarized fermionic atoms at zero temperature. The atoms 
interact via a short-range (tunable) attractive interaction. We demonstrate analytically a curious 
property of the gas at unitarity. Namely, the correlation energy of the gas, evaluated by second 
order perturbation theory, has the same density dependence as the first order exchange energy, 
and the two almost exactly cancel each other at Feshbach resonance irrespective of the shape of 
the potential, provided (fir B ) ^ 1. Here (/i) _1 is the range of the two-body potential, and r s is 
defined through the number density, n = 3/(47rrf). The implications of this result for universality 
is discussed. 



I. INTRODUCTION 

Consider a dilute gas of N ^S> 1 neutral fermionic atoms 
(mass M) at T — interacting with a short-range attrac- 
tive potential. In general, the properties of the dilute gas 
are determined by the number density n, and the scatter- 
ing length a. The Hamiltonian of this TV-particle system 
reads 



r,; - I\ 



(1) 



i<] 



Not written explicitly here, there is also an external po- 
tential Woxt(r) that forces the N atoms to stay within 
a large box with volume Q [where v ex t(r) = 0]. The 
attractive interaction potential is assumed to have the 
2-parameter form 



v(r) = -v f([ir) 



(2) 



where vq > is the strength of the interaction, R — f 
is its range, and f(x) is a dimensionless function. 

In the true ground state of the Hamiltonian (1) the 
attractive atoms may form dimers or even clusters. We 
are, however, looking for a metastable state where there 
is a dilute gas of separated atoms with uniform density n, 
satisfying the condition (jur s ) ^> 1 where n = 77 = 4^73- ■ 
Even then, for a weak vq, there will be BCS-type pair- 
ing, followed by dimer formation as the strength of the 
interaction increases. This was predicted long back by 
Leggett [1], and has been observed experimentally [2]. 
For the density functional analysis of the uniform gas 
at Feshbach resonance, we shall disregard the BCS con- 
densed pairs in this paper. 

To study the effect of the attractive interaction v(r), we 
consider the corresponding atom-atom scattering prob- 
lem in the relative s-state. Separating the center-of-mass 
motion, we are left with the relative Hamiltonian 



Keeping the range of the potential small enough such 
that (fJ,r s ) ^> 1, the strength vq is adjusted such that the 
potential can support a single bound state at zero energy. 
This happens when the scattering length a — > 00, leaving 
no length scale from the interaction. Such a tuning of 
the interaction is possible experimentally, and gives rise 
to Feshbach resonance [3]. The scattering cross section 
in the given partial wave (s-wave in our case) reaches the 
unitary limit, and the gas is said to be at unitarity. It is 
then expected to display universal behavior [4] . Note that 
at Feshbach resonance, there is no length scale left other 
than the inverse of the Fermi wave number kp, where 
kp = (Z^n) 1 / 3 . The energy per particle, E/N, as a 
function of the density n, should therefore scale the same 
way as the noninteracting kinetic energy, ^h 2 kp/2M oc 
n 2 / 3 . There has been much interest amongst theorists to 
calculate the properties of the gas in the unitary regime 
(kp \a\ 1). In particular, at T = 0, the energy per 
particle of the gas is calculated to be 



E_ 3 h 2 k 2 F 
N ~ 5 2M 



(4) 



where £ ~ 0.44 [5]. The experimental value of £ is about 
0.5, but with large error bars [6]. Recently, there have 
been two Monte Carlo (MC) finite temperature calcula- 
tions [7, 8] of an untrapped gas at unitarity, where various 
thermodynamic properties as a function of temperature 
have been computed. It is clear that at unitarity, the 
kinetic and potential energies should scale the same way. 
This has been assumed a priori in a previous density 
functional treatment of a unitary gas [9]. However, such 
a scaling behavior is not evident from the density func- 
tionals for the direct, exchange and correlation energies 
[10] (see sects. II and III). The aim of the present paper 
is to examine this point in some detail. In particular, we 
are able to show analytically that the leading contribu- 
tion of the correlation energy (calculated in second order 
perturbation theory), cancels the first order exchange en- 
rgy almost exactly at Feshbach resonance. This happens 
irrespective of the shape of the potential as specified by 
/(/zr), provided the condition (fir s ) ^> 1. We show that 
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our general Eq.(24) (derived later in the text) that en- 
sures such a cancellation is satisfied at unitarity for a 
variety of 2-parameter potentials, including the square 
well and the delta-shell, as well as the smoothly varying 
cosh~ 2 (/xr) and Gaussian potentials. This is the main re- 
sult of the present work. The implications of this result 
for universality is marginal. This is because these po- 
tential energy terms, in the limit of (fir s ) ^> 1, are very 
small compared to the kinetic energy [4]. For a moder- 
ately large value like (/zr s ) — 3 howevr, these terms are 
comparable in magnitude to the kinetic energy (sect. IV). 
Even then, the cancellation of the first order exchange, 
and the second order perturbativc terms leave the di- 
rect first order term in tact. In the electron gas, this 
(repulsive) term got cancelled by the interaction of the 
electrons with the positive ionic background. There is no 
such mechanism of cancellation here, unless we assume, 
rather arbitrarily, that the short-range interatomic repul- 
sion cancels this direct (attractive) contribution. Even 
without any such assumptions, however, our main result 
(Table I), applicable at Feshbach reonance, is interesting 
from the angle of potential theory. 



A. First order 

Formally, the first-order correction, 

=<Wnt|£0>, (7) 

has a direct contribution U{r s ,fi) — Nu(r s ,n) with 



u(r s ,») = f- d 3 r dV«(|r-r'|) 
3 v 



2 (fir s y 



(8) 



Here, f 2 = / °° dxx 2 f(x). 

The other first-order contribution is the exchange en- 
ergy E x (r s ,fi) = Ne x (r s ,fj,), [4] 

e x (r s ,^i) = drji{k F r) 2 v{r). (9) 

Here, ji(z) is a spherical Bessel function. Since v(r) is 
short-range and k F is small in a dilute gas, we can use 
the small- z expansion ji{z) = § + 0(z 3 ) to find 



II. PERTURBATION EXPANSION 



e x (r s ,M) = ^p^/ 2 + 0(^) 5 - 



(10) 



Treating the interaction (2) as a weak perturbation in 
the Hamiltonian (1), the unperturbed energy is the 
kinetic energy of a non-interacting Fermi gas, 



= Nt s (r s ) 



AT 3h 2 k 2 F 

N -. 

5 2M 



(5) 



Here, k F 



and 



* 3 = ± 

1 9ir 



. The corresponding ground 



state | #o) is a Slater determinant of plane waves. 

In terms of dimensionless coordinates Xj = /ir^, the 
Hamiltonian (1) can be written as 



M 

h 2 ^ 



H 



1 N 

= -^E v ?- a E/(n- x ^i)' a = 

i— 1 i<j 



Mv 



(6) 

This suggests that the perturbation parameter is not re- 
ally small at unitarity. For example, for the square-well 
potential, the zero-energy single bound state occurs when 
A = 2j- (see sect. III). Nevertheless the low-order terms 
can point to important information, even when the ex- 
pansion is divergent [11]. In our problem, there are three 
parameters, ji, vq, and r s . The unitarity condition re- 
lates jU and wo, so two independent parameters are left. 
One of these may be taken to be the small parameter 
C = (/ir s ) _1 . The remaining free parameter vq may be 
chosen independently of £ to fulfill the unitarity condi- 
tion. 



B. Second order 

1. General expressions 
Also the second-order correction, 



Em = _j- K*At|*b)| 2 = , 



n#0 



E n — Eq 



p ( 2 ) + p (2) 
°dir ' °ex 



, (11) 



has a direct and an exchange contribution [12], 



r d 3 fci d 3 k 2 

Jd q • (q + ki - 



k 2 )' 



64tt 5 \H 2 fj, 2 ) > 

x / d 3 h d 3 k 2 —. — / 

Jd q ■ (q + ki - k 2 ) 



(12) 



4 2 > s ,M)= + ^rs,o4/A/fe)x 



(13) 



While v 2 (2M/h 2 n 2 ) has the dimension energy, the inte- 
gration variables q, ki, and k 2 are dimensionless here. 
The domain of the integral over d 3 k\ d 3 k 2 depends on q, 

D: |k 1 |,|k 2 |<l; |ki+q|,|k 2 -q|>l. (14) 



3 



Furthermore, f(y) is a dimcnsionless transform of f(x), 



f(y) = / dxx f(x)j {yx) 
Jo 

1 f°° 

= - dxxf(x) sin(j/x). 

y Jo 



(15) 



To recover Eqs. (8) and (9) of Ref. [12], put M = m e , 
v = -e 2 fj,, and f(x) = \ or f(y) = ±r, such that v(r) 

' becomes the electronic Coulomb repulsion. (Note that 
Ref. [12] uses Rydberg units, m e e i /2h 2 = e 2 /2a B = 1.) 



2. The limit fir s > 1 

For a dilute gas (small kp) with short-range interaction 
(large /u), Eqs. (12) and (13) can be evaluated in the limit 
[i/kp = afj,r s S> 1 where a 3 — Following Ref. [13], 
we choose a number q\ such that 1 <C qi <C [ijkp and 
split the integrals over d 3 q into two parts, 



[ d 3 q = f d 3 q + f d 3 q. 
J J q<qi Jq>qi 



(16) 



In the first part with q < q\, we have q <C n/kp and 
|q + ki - k 2 | < fi/k F (note that |ki|, |k 2 | < 1 < qi). 
Therefore, we may expand f(y) = f2 + 0(y 2 ) in Eqs. (12) 
and (13) and keep the leading term / 2 only. The sum of 
the two resulting q < qi contributions reads 



e q<qi ( r S' A 4 ) — 

d 3 q 

q<qi 



2M 



64tt 5 \h 2 n 2 

rf 3 fcl d 3 k 2 

D q (q + ki - k 2 )' 



hi 

2«F ,2 
J ,A J2 



(17) 



The number q\ can be chosen independently of /i/kp 3> 
1, despite the condition 1 <C q\ <C [ijkp. Then, the 
integral in Eq. (17) is a finite constant and we conclude 
[13] 



(2 



: 9l (>s>M) = 0(pr s ) 4 . 



(18) 



In the second part q > qi » 1 of the integral (16), we can 
put q + ki — k 2 s=s q, since |ki|, |k 2 | < 1. The resulting 
contributions to Eqs. (12) and (13) add up to 



o (2) 
°q>qi 



3 /4tt 

'64tt 5 It 



/ 2M 



/r 4 



x/ ^/(^V (19) 
where / D d 3 /ci (i 3 fc 2 = (^) 2 has been used. Now, 

) 

*/(j) 2 (20) 



9><J1 



r~ 47r 



where yi = k F qi/fj, <C 1. If dyf(y) 2 in Eq. (20) did 
not depend on y\, expression (19) did rigorously have the 



order 0(pr s ) 3 . However, using the small-y expansion 
f(y) = h + 0(y 2 ), we have ^ dyf(y) 2 = f 2 Vl + 0(y 3 ). 
Consequently, shifting the lower limit y\ of the integral 
(20) to zero does not affect the leading-order contribution 
to expression (19), 



4>9i( r «'^) = °(M^) i - 



(21) 



Therefore, the quantity (18) does not contribute to the 
leading order of = + e$ which is purely due to 
expression (19), 



3 / 2M 



Air \ft 2 /i 2 / (/ir s ) 3 



U -°-F + 0(^r s )- 4 (22) 



where F = J °° dyf(y) 2 



III. DENSITY SCALING AT UNITARITY 

If the perturbation expansion is convergent [12], the 
total energy E(r S7 p) = Ne{r s ,p) of the gas can be ex- 
pressed in the form 

oo 

e(r 8 ,fj) = i s (r s ) + e x (r s ,^) + ^e( n) (r s ,/i). (23) 

At unitarity, when the relative Hamiltonian (3) has a 
single bound state at zero energy, the exchange plus cor- 
relation energy e x + YLn=2 e ^ should display the same 
density scaling as the kinetic energy, t s (r s ) cx r~ 2 cx p 2 / 3 . 
This is obviously not the case with any one of the present 
(leading-order) results (10) and (22). However, since the 
exchange energy (10) and the second-order correlation 
energy (22) have opposite signs, they can cancel each 
other at some value of fj,. This happens when 



Mv 



7t/2 

2F : 



(24) 



where / 2 = f£°dxx 2 f(x) and F = J^°dyf(y) 2 . This 
is the main result of our paper, and we check it by tak- 
ing four different potentials. The results of this analysis, 
summarized in Table I, are discussed in detail below. 

Generally, we need an eigenfunction ip(r) = ^p- 
of the relative Hamiltonian (3) with eigenvalue zero. 
Writing u(r) = (j>(pr), the corresponding dimcnsionless 
Schrodinger Equation reads 



4>"(x) = -\f(x)ct>(x), 



X 



Mv 



(25) 



Precisely, we wish to determine that particular value A u t y 
of A for which this zero-energy solution is the only bound 
state. Then, <p(x) must obey 0(0) = 0, <j)'{x) < for 
x > 0, and <f)(x) — > const, for x — > oo. In the following 
examples (A-D), the solution <f>{x) can be found analyti- 
cally or numerically. 
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(A) Square-well potential of radius Ro = 1/ fi: 
v(r) = -voO(Ro - r) , 



(26) 



where 0(z) denotes the Heavyside step function, 0(z) = 
1 for z > and 0{z) = for z < 0. By setting the 
dimensionless variable fir — x, we see that f(x) = 0(1 — 
x) . The square- well potential (26) supports a single zero 
energy bound state when the LHS of Eq.(24) is A u ty = 
7r 2 /4. It may be easily checked analytically that for the 
square- well potential (26), ft = | and F = so that the 
RHS of Eq.(24) is §, very close to its LHS, tt 2 /4 = 2.47. 

(B) Rosen- Morse hyperbolic potential [5] . This poten- 
tial is given by 



v(r) = —v sech 2 (fir) , 



(27) 



which suppotrs a single zero energy bound state when 
the LHS of Eq.(24) is A uty = 2 instead of tt 2 /4. For this 
potential, it is easy to check that = 7r 2 /12. The quan- 
tity F, however, has to be calculated numerically, and 
is given by F = 0.596. Again, Eq.(24) is approximately 
satisfied, since its RHS for this potential is 2.17. 

(C) Delta-shell potential [14]. Consider the potential 



v(r) 



Ro) 



-V W S(r 

——sf— 

^MRo \R 
-v f(fir) . 



(28) 



Thus, we have i>o 
So we get f 2 = 1, f(y) = 



'I MRo ' 



/' = 3^, and f(x) = 5(x-l). 

Hence the 



-Ro 

and F 



sin y 

y ' ^ — 2 ' 

RHS of Eq. (24) is unity. The LHS is (r)R ), which is 
exactly unity when the s-state scattering length goes to 
infinity [14]. Thus Eq.(24) is exactly obeyed in this case. 
(D) Gaussian Potential. 



v(r) = — «ocxp(— fi r ) 



(29) 



For this example, f(x) = cxp(— x 2 ) in Eq. (2). We find 
h = \ypH and F = |(f ) 3/2 so that the RHS of Eq. (24) 
becomes irf2/2F = 2 3 / 2 . Solving Eq. (25) numerically 
for this f(x), we obtain a single bound state at zero 
energy when the LHS of Eq. (24) is A uty = 0.949 x 2 3 / 2 , 
close to 2 3 / 2 . 

Note, however, that contributions 0(fir s )~ 3 may also 
come from higher order terms of the perturbation expan- 
sion in section II, since that expansion is carried out with 
respect to the parameter A = Mva/h 2 fi 2 7 but not 1 / fir s . 



IV. DISCUSSION 

The dimensionless Hamiltonian h from Eq. (6) depends 
on the dimensionless paramaters 



A = 



Mv 



(30) 



TABLE I: The moments fc and F of four different profiles 
f(x) for the potential (2). A u t y is the value at unitarity of the 
parameter A in Eq. (25). At unitarity, the ratio Q of the LHS 
of Eq. (24) to the RHS is always close to 1. 



/(*) 


h 


F 


Auty 


Q 


9(1 -x) 


i 

3 


15 


4 


0.987 


sech(x) 2 


vr 2 
12 


0.596 


2 


0.922 


5(1 -x) 


1 


7T 
2 


1 


1.000 


exp(—x 2 ) 




8^2 > 


2.684 


0.949 



and, not written explicitly, x s — fir s . The perturbation 
expansion of the ground-state energy of h reads 



e(x s , A) 



n=0 



e n (x s )y 



(31) 



The ground-state energy of the original Hamiltonian H, 
with three independent parameters, is then given by 



E(r s ,fi,X) 



h 2 fi' 

M 

*0 n OO 

h l fi z 

M 



e(fir s ,\) 

OO 

J2 £ n(^ir s )X n . 



(32) 



n=0 



For fir s 3> 1, we may expand 



e n (fJ-r s ) = 



r (vr s 



(33) 



From Eq. (5), we have £02 = 2 while e 0m = for 

m/2. Eqs. (8) and (10) imply that E\ m = for m < 3 
and ei3 = A(-§ + |)/ 2 . Eventually, due to Eq. (22), 
£ 2m = for m <C 3 and £ 23 = N(--^)2F. 

So far as the unitary point is concerned, we are inter- 
ested in a situation where kp\a\ >1> IzfRo ~ (firs)^ 1 . 
In view of the fact that the perturbation series above does 
not converge at unitarity, how significant is our low order 
perturbation calculation in this situation ? Note that our 
first order direct and exchange (potential) energy terms 
given by Eqs. (8,10) are the same as those obtained in 
the Hartree-Fock calculation (see, for example, Eq.(10) 
of Heiselberg [4] ) . How big are these terms at unitarity 
compared to the kinetic energy per particle ? Taking the 
example of the square-well potential discussed earlier, it 
is straight forward to show that our exchange term (10) 
at Feshbach resonance is 



, . tt /9tt\ 1/3 E f 
e * (rs ' M)= 18UJ J^s- 



(34) 



■5 



For the square- well example, 



(kpa) = 



9tt 



1/3 



1 - 



tan\/A 



(35) 



At unitarity, the RHS diverges for any finite value of 
(//r s ), how ever large. Even in the neighbourhood of 
unitarity, it is possible to have (kp\a\) S> 1 for (/Ur s ) S> 1. 
From Eq.(34), we note that too large a choice for fjur s ) 
would make e x negligible against §-Ef- Instead, taking 
a modestly large value, /ir s = 3, we obtain the ration 
of e x to kinetic energy per particle to be about 0.56. 



Noting that e x has a different density-dependence than 
the kinetic energy per particle, its cancellation with the 
second order perturbative correlation term helps towards 
scale invariance, but only if there is a mechanism for the 
direct first order term to be cancelled. 

We conclude by emphasizing that the new result in this 
paper is displayed in Table 1, and should be of interest 
from the point of view of potential theory. 

The authors would like to thank Brandon van Zyl for 
discussions. This research was financed by NSERC of 
Canada. 
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